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Abstract 

Giardia duodenalis is a common and widespread intestinal protozoan parasite of both humans and animals. Previous 
epidemiological and molecular studies have identified Giardia infections in different animals and humans, but only limited 
information is available about the occurrence and genotypes of Giardia in cattle in China. In this study, we determined the 
occurrence of giardiasis and genetically characterized G. duodenalis in dairy cattle in Henan Province, central China. The 
overall prevalence of G. duodenalis was 7.2% (128/1777) on microscopic analysis, with the highest infection rate (22.7%) in 
calves aged less than 1 month. G. duodenalis assemblages and subtypes were identified with multilocus genotyping based 
on the SSU rRNA, p-giardin (bg), glutamate dehydrogenase (gdh), and triosephosphate isomerase [tpl) genes. Two 
assemblages were detected in the successfully sequenced samples: assemblage A (n = 58), assemblage E (n = 21), with a 
mixed E and A assemblage (n = 2). Four novel subtypes of the gdh gene and seven of the bg gene were found among the G. 
duodenalis assemblage E isolates. Using the nomenclature for the multilocus genotype (MLG) model, nine novel multilocus 
genotypes E (MLGs E1-E9) and three MLGs A (a novel subtype Al, previously detected subtype AII-1, and a combination of 
both) were identified. MLG AII-1 identified in this study may be an important zoonotic subtype. The dairy cattle in Henan are 
a potential public health concern. 
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Introduction 

Giardia duodenalis (syn. G. lamblia, G. intestinalis) is a common 
protozoan in cattle worldwide. It often presents no or mild 
symptoms in adult cattle, but clinical manifestations of diarrhea, 
weight loss, and malabsorption can present in calves [1]. Previous 
epidemiological studies have shown that cattie are the main source 
of G. duodenalis infection in humans, contracted by the fecal-oral 
route or by the ingestion of contaminated food or water [2]. 
Therefore, the role of cattle as reservoirs of G. duodenalis and its 
potential threat to public health are of increasing concern [3] . 

The accurate identification of parasites is central to the effective 
control of parasitic diseases. Molecular studies have confirmed that 
G. duodenalis is a species complex, comprising eight distinct 
assemblages/genotypes (A-H) [4,5] that appear to have different 
host ranges. Of these assemblages, only A and B infect humans 
[6,7]. So far, the zoonotic assemblage A and B and the livestock- 
specific assemblage E have been detected in catde [8-16], 
Assemblage E is the predominant genotype in most countries, 
including Belgium, United States, Canada, Denmark, Australia, 
and Portugal [8- 1 3] . In contrast, assemblage A or B is occasionally 



reported to be the most common genotype in Italy, Canada, and 
New Zealand [14-16]. 

In China, the occurrence and molecular studies of Giardia have 
been reported in rabbits (subtypes B-I and B-VIII), monkeys 
(subtypes A-II and B), macaques (subtypes A-II and B), sheep and 
goats (genotypes A, B, and E), dogs (genotypes C, D, and A), and 
humans (subtype A-I, A-II, and B) [17-23]. However, there has 
been only one study of G. duodenalis infection in calves in 
Heilongjiang, China [24]. Therefore, because there are relatively 
few prevalence data available on G. duodenalis in cattle in China 
and an even greater lack of molecular data, we determined the 
prevalence of G. duodenalis in dairy cattle in Henan, China, and 
characterized it at the molecular level, using multilocus genotyping 
at the SSU rRNA, fi-giardin (bg), glutamate dehydrogenase (gdh), 
and triosephosphate isomerase (tpi) loci. 

Materials and Methods 

Ethics statement 

This study was performed strictly according to the recommen- 
dations of the Guide for the Care and Use of Laboratory Animals 
of the Ministry of Health, China. Our protocol was reviewed and 
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approved by the Research Ethics Committee of Henan Agricul- 
tural University. All fecal specimens were collected from animals 
with the permission of the farm owners. No specific permits were 
required by the authorities for specimen collection. The field 
studies did not involve endangered or protected species. 

Sample collection and microscopy 

From August 2008 to October 2009, 1777 fecal samples were 
randomly collected from dairy cattle aged from newborn to 2 years 
old on 15 different intensive farms in Henan Province, central 
China (Figure 1, Table 1). New fecal specimens from cattle at 
Farm 1 were also collected monthly for 1 year to investigate the 
seasonal variations in G. duodenalis infections. Farm 1 is the largest 
dairy farm in Henan Province, located in a suburb of the city of 
Zhengzhou, consisting of approximately 1000 animals, including 
calves (<6 months old), heifers (6-24 months old), and cows (>2 
years old). The farm ranked among the top producing dairy farms 
in Henan. All dairy cattle were housed in different age groups in 
free stalls indoors. 

The fresh feces were placed into clean plastic bags marked with 
the date, age, and geographic origin, transported immediately to 
the laboratory, stored at 4°C, and then examined with a Lugol- 
staining technique and microscopy at 400 x magnification. All G. 
duodenalis-positive samples were stored in 2.5% (w/v) potassium 
dichromate solution at 4°C for molecular characterization. 

DNA extraction and PCR amplification 

The samples were washed three times in distilled water with 
centrifugation at 3000 x g for 10 min to remove the potassium 
dichromate. The genomic DNA was extracted from the G. 
duodenalis samples using the E.Z.N.A. Stool DNA kit (Omega 
Biotek Inc., Norcross, GA, USA), according to the manufacturer's 
instructions. The extracted DNA was eluted in 100 |il of AE 
elution buffer and stored at — 20°C. 

The genotypes of the G. duodenalis samples were determined with 
nested PCR amplification of the SSU rRNA, bg, gdh, and tpi genes, 
according to previous studies [10,25-27], with several modifica- 
tions. The primers used in the PCR analysis of all gene targets, the 




Figure 1. Specific locations at which samples were collected in 
this study. ▲ locations. 

doi:1 0.1 371 /journal.pone.01 00453.g001 



annealing temperatures, and the sizes of the expected PCR 
products are listed in Table 2. The PCR reactions for the bg, gdh, 
and tpi loci were conducted in 25 |iL reaction mixtures containing 
of 1 x PCR buffer (TaKaRa Shuzo Co., Ltd., Otsu, Japan), 
200 |j,M each dNTP (TaKaRa Shuzo Co., Ltd.), 0.4 uM each 
primer, 1 unit of TaKaRa rTaq DNA polymerase (TaKaRa 
Shuzo Co., Ltd.), and 2 uL of DNA sample. In the SSU rRNA 
protocol, 1 x GC buffer II (TaKaRa Shuzo Co., Ltd.) and LA Taq 
DNA polymerase (TaKaRa Shuzo Co., Ltd.) were used instead of 

1 x PCR buffer and rTaq. For the second amplification of bg, 
10 uL of the first PCR amplicon was diluted in 90 |0,L of water and 

2 uL of that dilution was used as the template for the second 
amplification. DNA samples from G. duodenalis isolates from 
humans and distilled water were used as the controls for each 
target-gene-based PCR analysis. The PCR products were visual- 
ized on a UV transilluminator after electrophoresis in 1.5% 
agarose gels and staining with ethidium bromide. 

Sequence analysis 

All nested-PCR amplicons were sent to Beijing Nuosai 
Biological Engineering Biotechnology Company for two-direc- 
tional sequencing on an ABI PRISM 3730 XL DNA Analyzer 
(Applied Biosystems, USA). The G. duodenalis genotypes were 
identified by alignment with reference sequences downloaded 
from GenBank (http://www.ncbi.nlm.nih.gov) using MEGA 4 
[28,29]. A phylogenetic analysis was performed using the 
neighbor-joining and maximum composite likelihood methods 
carried out with Tamura-Nei model. The consensus tree was 
constructed after bootstrap analysis with 1000 replications. 

Statistical analysis 

The % 2 test was used to compare the prevalence of G. duodenalis 
between different cities, the association between infection rate, 
age, and sampling season, and between the assemblage distribu- 
tions and age. Differences were considered significant at P<0.01. 

Nucleotide sequence accession numbers 

All nucleotide sequences were submitted to the National Center 
for Biotechnology Information (NCBI) GenBank database under 
the following accession numbers: KF843921-KF843922 for the 
SSU rRNA gene, KF843923-KF84393 1 for the gdh gene, 
KF843932-KF843940 for the bg gene, and KF843941- 
KF843948 for the tpi gene. 

Results 

Prevalence of G. duodenalis 

The total prevalence of G. duodenalis was 7.2% (128/1777); 13 of 
15 farms were Giardia-positive and the infection rate (95% CI) 
ranged from 0% to 14.4% on different farms (Table 1). The 
highest infection rate was in Luoyang (11.6%, 28/241), followed 
by Zhengzhou (7.2%, 72/1002), Shangqiu (5.8%, 9/154), 
Xinxiang (4.9%, 9/184), and Jiaozuo(4.1%, 10/244). Among the 
different age groups of cattle, calves aged less than 1 month had 
the highest infection rate (% = 171.1, P<0.01), whereas no Giardia 
cysts were detected in adult cattle older than 1 year (Table 1). 

The association between the sampling season and the preva- 
lence of G. duodenalis was assessed using fecal samples from Farm 
Zz-1 (Figure 2). The prevalence was highest in February and 
lowest in June. The animals were at significantly greater risk of 
Giardia infection in winter (November-February) (% 2 = 29.34, P< 
0.01). 
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Table 2. Target, primers, amplicon size, annealing temperature, and 


main use of the four G. duodenalis genotyping loci. 




gene Primer (sequence 5 -3 ) 


Fragment 
length (bp) 


Annealing 
temperature CO 


Usage(s) 


References 


16S rRNA Gia2029 (AAGTGTGGTGCAGACGGACTC) 


292 


55 


genotyping 


9 


Gia21 50c(CTGCTGCCGTCCTTGGATGT) 


RH 11 (CATCCGGTCGATCCTGCC) 




59 






RH4(AGTCGAACCCTGATTCTCCGCCCAGG) 


bg G7(AAGCCCGACGACCTCACCCGCAGTGC) 


384 


65 


Genotyping and 
subtyping 


24 


G759(GAGGCCGCCCTGGATCTTCGAGACGAC) 


G759(GAGGCCGCCCTGGATCTTCGAGACGAC) 




65 






G376(CATAACGACGCCATCGCGGCTCTCAGGAA) 


gdh Ghdl (TTCCGTRTYCAGTACAACTC) 


520 


50 


Genotyping and 
subtyping 


25 


Gdh2 (ACCTCGTTCTGRGTGGCGCA) 


Gdh3 (ATGACYGAGCTYCAGAGGCACGT) 




50 






Gdh4 (GTGGCGCARGGCATGATGCA) 


tpi AL3543 (AAATIATGCCTGCTCGTCG) 


530 


50 


Genotyping and 
subtyping 


26 


AL3546 (CAAACCTTITCCGCAAACC) 


AL3544 (CCCTTCATCGGIGGTAACTT) 




50 






AL3545 (GTGGCCACCACICCCGTGCC) 



doi:1 0.1 371 /journal.pone.01 00453.t002 



Molecular identification and polymorphisms of G. 
duodenalis isolates 

We genotyped 128 G. duodenalis-positive samples based on four 
loci, and 73 SSU rRNA, 67 tpi, 11 bg, and 61 gdh gene sequences 
were obtained. Fifty-six samples were successfully sequenced at all 
four genes, whereas the remaining samples were amplified at only 
one to three genes (Table 3). All nucleotide sequences obtained in 
the present study were aligned with G. duodenalis reference 
nucleotide sequences from GenBank and analyzed using ClustalX 
1.83. Two G. duodenalis assemblages were identified: assemblage E 
(n = 58), reported to infect only hoofed livestock, and assemblage A 
(n = 21), which infects humans and a number of other mammals. 

The genetic diversity of G. duodenalis among assemblage E and A 
was observed at the tpi, bg, and gdh loci. Six different subtypes of 



assemblage E were identified at the tpi locus (KF843941- 
KF843946), and these showed 100% similarity to sequences 
available in GenBank with accession numbers of AY655705, 
EF654683, JN162351, AB569406, EF654689, and JF792419, 
respectively. Seven subtypes of assemblage E were also identified 
at the bg and gdh loci. Four subtype sequences at the gdh locus 
(KF843926-KF843929) have not been reported previously, and 
showed similarity of 99% to assemblage E (accession 
no. EF507645, from cattle in Brazil). The remaining three subtype 
sequences at the gdh locus (KF843923-KF843925) showed 100% 
similarity to isolates B7 from cattle (EF507645), SG-10 from sheep 
(KC960647), and Ca39 from cattle (AB569388), respectively. 
Subtypes AI and All of assembly A were identified at the tpi, bg, 
and gdh loci, with subtypes AI and All identical to sequences with 




Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec 



Month 



Figure 2. Prevalence of G. duodenalis infections in dairy cattle in different months. 

doi:1 0.1 371 /journal.pone.01 00453.g002 
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Table 3. Genotype distributions of the 16s rRNA, bg, gdh, and tpi gene sequences. 





SSU rRNA 


tpi 


bg 


gdh 


Total 


E 


E 


E 


E 


38 




NR 


E 


NR 


3 




NR 


NR 


E 


3 




NR 


E 


NR 


8 


NR 


E 


E 


NR 


6 




A 


E 


A 


2 


A 


A 


A 


A 


18 




A 


NR 


NR 


1 




A 


A 


NR 


2 



NR: no result. 
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accession numbers AF069556 (subtype AI) and AF069557 
(subtype All) at the tpi locus, AY655702 (subtype AI) and 
AY072723 (subtype All) at the bg locus, and EF507642 (subtype 
AI) and EF507674 (subtype All) at the gdh locus. The variations 
among the new subtypes at different nucleotide sites are 
summarized in Table 4. 

Distribution of G. duodenalis assemblages 

Across all farms, 72% of the G. duodenalis-positive samples were 
infected with assemblage A, while 26% were infected with 
assemblage E. Assemblage E was found on twelve of fifteen farms, 
while assemblage A was only found on three farms: two of three 
farms being infected with assemblage E and assemblage A, one 
being present as a single assemblage A on Xx- 1 (The data was not 
shown). Fifty-eight G. duodenalis assemblage E and twenty-one G. 
duodenalis assemblage A isolates were assigned to different ages 
groups. Assemblage E was found in all age groups and the highest 
percentage was in 61-90-day-old cattle (% 2 = 25.02, P<0.01), 
while assemblage A was only observed in calves < 1 80 days old, 
with the highest percentage in calves <30 days old (% 2 = 29.58, P< 
0.01; Figure 3). 

Multilocus genotyping 

Because there was no variability in the SSU rRNA gene among 
assemblage A or assemblage E isolates, the bg, gdh, and tpi loci were 
used to determine the G. duodenalis subtypes according to the 
established nomenclature, based on multilocus sequence polymor- 
phisms. Fifty-six isolates were successfully subtyped at all three 
genes, forming nine different assemblage E MLGs and three 
assemblage A MLGs (Table 4). Of the assemblage E MLGs, the 
greatest number of isolates belonged to MLG El (n = 16), followed 
by MLG E2 (n = 8). The other seven assemblage E MLGs were 
represented by one to five isolates. Of the assemblage A MLGs, 
five isolates belonged to subtype All (a known MLG AII-1), 12 
isolates were subtype AI (a novel MLG AI), and only one isolate 
showed a combined MLG AII-1 and novel MLG AI subtype. The 
GenBank accession numbers of the reference MLGs (AI-1, AI-2, 
AII-1, AII-2, and ASSS) and all the MLGs identified in this study 
are listed in Table 5. To clarify the genetic relationships between 
the different MLGs, a phylogenetic analysis was performed based 
on a concatenated dataset of bg, gdh, and tpi sequences. The 
different MLGs of G. duodenalis were included as in-groups for 
improved resolution and a better topology of the evolutionary tree. 
All MLGs in the current study were assigned as assemblage A and 
assemblage E. In the assemblage A MLGs, one subtype clustered 



with the reference MLG All- 1 , whereas the novel subtype AI and 
the combined MLG clustered with the reference MLG AI- 1 . The 
MLGs E1-E9 clustered broadly with previous assemblage E 
MLGs (Figure 4). 

Discussion 

In this study, the overall infection rate of G. duodenalis in cattle, 
ranging from newborn to 2 years old, was 7.2% (128/1777), which 
is a little higher than that recently reported in same-aged cattle in 
Heilongjiang Province, China (6.4%, 41/643) [24]. As is well 
known, the determinants of infection rates are complex and the 
rates are often affected by many factors, including the ages of the 
animals, sample size, examination methods, different management 
systems, the timing of specimen collection, and geoecological 
conditions. Therefore, it is difficult to explain the actual difference 
between the prevalence of G. duodenalis in this study and those in 
other countries, such as the United States (44%), Canada (42.0%), 
Galicia (30.1%), Denmark (43.6%), Belgium (31.3%), Vietnam 
(10.2%), and New Zealand (31.0%) [8,11,15,30-33]. 

There was a significant association between the age of sampling 
and the likelihood of infection. Results on the prevalence of G. 
duodenalis in catde in the current study revealed that infection was 
frequent in calves younger than 6 months old, but was rare in 
adult cattle older than 1 year. This finding was in accordance with 
previous reports regarding giardiasis as a common infection in 
immature farm animals [24,34,35,36,37]. The highest prevalence 
was observed in calves less than 1 month old (22.7%), which was 
different from previous reports. Huetink et al. reported that G. 
duodenalis cysts were rare in calves younger than 1 month and 
found the highest prevalence in 4-5-month-old animals [38]; 
Mark-Carew et al. [37] reported that young stock within the age 
range of 31-60 days were at the highest risk; The highest infection 
rate was in calves aged 2-3 months in O'Handley et al study [39]. 
Due to the lack of epidemiological data on bovine giardiasis, we 
could not find the true reasons for the relatively low age at which 
cysts had the highest prevalence in the current study. It might be a 
consequence of the calves management, the health status such as 
diarrhea or not, with or without the presence of other parasites 
such as Cryptosporidium, and so on. 

In this study, the peak G. duodenalis infection was noted in winter, 
which is identical to a previous study conducted in the Netherlands 
[38]. In that study, Huetink et al. repeatedly sampled a dairy farm 
over 1 year and found that the prevalence of G. duodenalis shedding 
peaked in December and February. They explained that in this 
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■ Assemblage E 
□ Assemblage A 




<30day 31-60 day 61-90 day 91-180 day 

Age group 



>180 day 



Figure 3. Frequency distribution of G. duodenalis assemblage E and assemblage A across different age groups. Sixteen, thirteen, 
twenty-two, six, and one of fifty-egiht G. duodenalis assemblage E isolates were grouped into calves aged <30 days, 31-60 days, 61-90 days, 91-180 
days, and >1 80 days, respectively. Twelve, four, two, three, and zero of twenty-one G. duodenalis assemblage A isolates were grouped into the 
corresponding groups. 
doi:1 0.1 371 /journal.pone.01 00453.g003 



period, because the weather is cold and the animals are housed 
indoors, the risk of infection with serious respiratory illnesses 
increases, making the cattle more susceptible to all kinds of 
infections, including G. duodenalis, as their immunity declined. 
These infected animals excreted G. duodenalis oocysts and then 
possibly infected other cattle in the cattle sheds via another vector: 
caretakers, cats, or mice. Similar climates might explain the 
consistent observations between the Netherlands and this study. 

Sequence analyses indicated that most of the dairy cattie were 
infected with livestock-specific G. duodenalis assemblage E (72%), 
whereas a minority harbored the zoonotic assemblage A (26%). 
These findings are consistent with previous reports from China, 
Belgium, United States, Canada, Denmark, Australia, and 
Portugal [8-13,24]. Interestingly, it was difficult to assign two 
isolates to an assemblage unequivocally because the genotyping 
results at the bg gene (assemblage E) were not consistent with those 
at the other two genes (assemblage A) (Table 3). This may indicate 



the presence of mixed genotypes. Only a few studies have reported 
mixed infections of G. duodenalis in cattle and other animals 
[7,26,40,41]. This discrepancy at different loci has important 
implications for molecular epidemiological studies, and using only 
one maker to assign an isolate to a specific assemblage is not 
always reliable. In contrast, multilocus genotyping is a suitable tool 
with which to interpret the molecular epidemiology of G. duodenalis 
infections. 

In a study of Cryptosporidium species and genotypes, the zoonotic 
species, C. parvum primarily infected calves under 2 months of age, 
while the species and genotypes that were not infectious for 
humans primarily infected post-weaned calves [42]. To identify 
whether this change occurs with G. duodenalis, we observed the 
differences in the assemblage distributions in the various age 
groups. The potentially zoonotic assemblage A was only found in 
calves < 1 80 days of age, whereas assemblage E was found in all 
age groups. The highest infection rates of assemblages A and E 
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Figure 4. Phylogenetic relationships of G. duodenalis MLGs. The phylogenetic tree was constructed using a concatenated dataset of the bg, 
tpi, and gdh gene sequences, and a maximum likelihood analysis and neighbor-joining analysis generated identical topologies. Sequences from this 
and previous studies are included in the analysis. Bootstrap values >50% are shown. G. duodenalis MLGs identified in this study are indicated by black 
triangles. 
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Table 5. Multilocus characterization of G. duodenalis isolates from dairy cattle based on the sequences of the bg, gdh, and tpi 
genes. 





Isolates (numbers) 


Reference sequences for bg, gdh and tpi genes 


MLG types 




Genotypes 


GenBank accession nos. 




Fy32 a (5) 


All All All 
All, Mil, Mil 


AVnV^VJ^ PP^n7fi7Zl 1 K7RQ7 
AYU/Z/ZD, LrDU/O/^f, Uj/o7/ 


All-1 b 


Xx1 a {12) 


Al Al Al 


iVfi'i^m AR1^Q7Q^ imnn 

A I ODD / UZ, AD I Dy /yD, LUZ I ZU 


A novel b 


Fy16 (1) 


All Al Al 
All, Al, Al 


Avn797'?^ AR1^Q7Q^ imnn 

AYU/Z/ZD, AD I Dy /yD, LUZ I zu 


Mixed" 


Swesheep055 


Al Al Al 


Al^floD, LrDU/OlU, LUZ I ZU 


AI-1 


ISSGCat4 


Al Al Al 


nDtUJJOJ, IVI0 1 +OU 1 +, nDJUJJOJ 


AI-2 


ISSGd168 


All All All 
All, All, All 


AVfl7'?7'?'5 PP^n7ft7A I K7RQ7 
AlU/Z/ZD, LrDU/O/^r, UD/Oy/ 


AII-1 


ISSGd107 


All All All 
All, All, All 


AVfl7'?7'?A PI n7R^nR I K7RQ7 
AlU/Z/Z'f, LUZ/oDUo, UD/Oy/ 


AII-2 


Swecat171 


AIM AIM AIM 
AIM, AIM, AIM 


nOfi^flfiAQ PI lfi^7c;RT PII7Rin(T? 
UL^DDUD'fy, LUDD/DOZ, LU/olUUZ 


AIII-1 


Swesheep069 


E E E 


nnnfifiiA nnifi^fin^ Pii7RiniQ 

L*M I IDDZ'f, UVJIOZDUD, LU/OIUiy 


MLGEI 


Swesheep026 


E,E,E 


EU769215, DQ182605, EU781019 


MLGEII 


Ly7 a (16) 


E3,E3,E1 


KF843934, KF843925, KF843941 


MLGEI b 


Ly10 a (8) 


E1,E3,E1 


KF843932, KF843925, KF843941 


MLGE2 b 


Zm68 a (5) 


E2,E1,E1 


KF843933, KF843923, KF843941 


MLGE3 b 


Sq51 a (2) 


E3,E2,E1 


KF843934, KF843924, KF843941 


MLGE4 b 


Ly24 a (2) 


E3,E1,E3 


KF843934, KF843923, KF843943 


MLGE5 b 


Dy61 a (2) 


E5,E1,E1 


KF843936, KF843923, KF843941 


MLGE6 b 


Ly34 (1) 


E4,E2,E1 


KF843935, KF843924, KF843941 


MLGE7 b 


Ly13 (1) 


E6,E3,E6 


KF843937, KF843925, KF843946 


MLGE8 b 


Jz37 (1) 


E1,E3,E4 


KF843932, KF843925, KF843944 


MLGE9 b 



a : Pooled sample. 
b : Present study. 

doi:1 0.1 371 /journal.pone.01 00453.t005 



were found in animals of 1 and 3 months old, respectively. These 
results are similar to those of previous reports in which the 
prevalence of assemblage E was highest in post-weaned calves and 
assemblage A was only detected at 6 months of age [9,30]. 
However, our findings differed from those of another study, in 
which assemblage E peaked in preweaned calves and assemblage 
A was detected in post-weaned calves and heifers [43] . 

The MLG model was developed by Caccio et al. in 2008 and 
used later by Sprong et al. and Lebbad et al. [7,26,44,45]. Using 
MLGs, we can better characterize the G. duodenalis in humans and 
animals from different geographic regions. In this study, three 
MLGs were identified in 18 assemblage A isolates, including a 
known MLG AII-1 (n = 5), a novel MLG Al (n=12), and a 
combination of both (n = 1). The distribution of the G. duodenalis 
multilocus subtypes is consistent with previous findings that 
subtype Al has a preference for livestock, whereas subtype All is 
commonly detected in humans [26,46]. The known MLG subtype 
AII-1 identified in this study is identical to a cat-derived strain 
isolated in Italy, and human-derived isolates in Italy, Sweden, and 



China [26,45,47]. The observation of a genetically identical G. 
duodenalis subtype in humans and animals in China and other 
countries indicates that MLG All- 1 may be an important zoonotic 
multilocus genotype. 

In conclusion, this study has reported the prevalence of Giardia 
in dairy cattle in Henan Province, China, and provides some 
preliminary data on the genetic diversity of the parasite in this 
region using multilocus genotyping. Whether MLG AII-1 is the 
main zoonotic subtype must be investigated further with 
systematic molecular epidemiological analyses of both humans 
and animals. 
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